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BalloonAscent is a new user-friendly software tool, which simulates the vertical and 
horizontal motions of high-altitude balloons such as are flown by the National Scientific 
Balloon Facility. The formulations are generalized as to allow simulated balloon flight in 
other atmospheres on different planets. This paper will describe the theoretical basis and 
how the tool is used to predict ascent behavior and float behavior under varying 
environmental conditions. 
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Nomenclature 

mass of lifting gas, kg 
temperature of lifting gas, °K 
pressure of lifting gas, Pascals 

specific gas constant of the lifting gas (2078.5 J/kg-°K for helium) 

specific heat at constant volume of the lifting gas (3 1 2 1 .5 J/kg-°K for helium) 

specific heat of the film material, J/kg-°K 

mass of film material in the inflated bubble, kg 

gas density, kg/m 3 

ambient air density, kg/m 3 

ambient air temperature, °K 

ambient air pressure, Pascals 

surface air pressure, Pascals 

air mass correction factor for smoke or haze in the atmosphere 
solar elevation angle above the true horizon, radians 
acceleration of gravity, m/s 2 

= design maximum volume of the gas envelope, m 3 
design maximum gore length, m 
orbital eccentricity 

mean orbital radius, astronomical units (earth R AU = 1.0) 

balloon skin bulk temperature, °K 

Stephan-Boltzman constant, 5.67E-8 Watts/m 2 -°K 4 

bulk averaged balloon skin absorptivity, visible sunlight 

bulk averaged balloon skin absorptivity, IR 

bulk averaged balloon skin IR emissivity 

bulk averaged balloon skin transmissivity, visible sunlight 

bulk averaged balloon skin transmissivity, IR 

bulk averaged balloon skin reflectivity 

direct solar flux, Watts/m 2 

diffuse planetary albedo flux, Watts/m 2 

diffuse up-welling planetary infrared flux, Watts/m 2 

diffuse down-welling infrared flux from the sky. Watts/m 2 
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I. Introduction 

T HE BalloonAscent balloon flight simulation code represents a ffom-scratch development using Visual Basic 5 
as the software platform. The simulation code is a transient analysis of balloon flight, predicting the skin and 
gas temperatures along with the 3-D position and velocity in a time and spatially varying environment. There are 
manual and automated controls for gas valving and the dropping of ballast. Also, there are many handy calculators, 
such as appropriate free lift, and steady-state thermal solutions with temperature gradients. 

The strength of this simulation model over others in the past is that the infrared environment is deterministic 
rather than guessed at. The ground temperature is specified along with the emissivity, which creates a ground level 
IR environment that is then partially absorbed as it travels upward through the atmosphere to the altitude of the 
balloon. 

The environment is specified with: 

a) Five-point Atmosphere model; temperature as a function of altitude which implies pressure 

b) Winds aloft; wind speed and direction as a function of altitude 

c) Launch location and month/day/time of launch 

d) Day/night ground temperatures and emissivity as a function of time 

e) Cloud altitude implying cloud temperature 

f) Ground albedo as a function of time 

g) Cloud albedo as a function of time 

h) Cloud sky obscuration fraction a function of time 

i) Sun elevation angle as a function of local time and latitude/longitude 

The high-altitude scientific balloon is a thermal vehicle. It lives in a mostly radiant thermal balance from the many 
heat sources that influence it. When warm, the gas expands and either pressurizes the gas envelope or is vented. If 
vented as happens with a zero-pressure balloon, then nighttime cooling makes ballast drops necessary. Soon one 
runs out of both ballast and gas, which is why performance predictions are critical. Super-pressure balloon design 
also benefits, as the day/night pressure drop is a critical bit of information during the design process. 

These are the basic heat sources that influence balloon flight: 

a) Direct visible light sunshine on the skin membrane 

b) Reflected diffuse sunshine in the form of planetary ground albedo 

c) Reflected diffuse sunshine in the form of planetary cloud albedo 

d) Diffuse infrared from the ground/atmosphere 

e) Diffuse infrared from the clouds 

f) Internal infrared “self-glow” of interior skin 

g) Atmospheric convection 

h) Diffuse infrared from the sky (important only in environments like Venus) 

To give one an idea of the influences at 33 kilometers altitude on an average sunny day, the fraction of heat loads 
absorbed in the skin of a large scientific balloon (10 to 40 million cubic foot volume) is approximately: 

Planetary surface IR 49% 

Direct solar 35% 

Indirect solar (albedo) 20% 

Self-glow internal IR 3% 

Atmospheric convection -7% 


II. Theoretical Development 

I will not describe every equation in the program, but rather focus on some of the more unique or essential 
formulations. 

A. Approximate Geometric Properties 

Depending on the three classes of balloons treated in the software (spherical, zero pressure natural shape, or 
super-pressure pumpkin), geometric properties are approximated with the following equations. 
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First, the unconstrained volume of the bubble, m 3 : 

Volume = M 


R„ 


gas gas 



Equations for a zero-pressure shape, closest match for a sigma shape = 0.2: 

Balloon diameter, top view, m: 


Diameter = 2.23 • 


0-75 mV 

v 71 J 


Length of gore exposed in the bubble, m: 


LgoreB = 1 -37 • Diameter 

Surface area showing in a zero pressure shape: 

A sur f = 2.582 ■ Diameter 2 

or 2 

A S urf = 4-94 • Volume * 


( 1 ) 



(4a, 4b) 


Maximum “crenellated” surface area in a zero pressure shape not fully inflated: 

2 

Aurf I = 4 -94 • Volume design 3 • (l - cos(^ • L goreB / L goreDesign )) (5) 


To include the effect of the crenellated excess surface area for convection and film mass calculations, an estimate is 
made regarding the effective exposed surface area: 

A e jfective = * A surf + 0.35 * A su ^ x (6) 

Top projected area: 

Ao P =^' Diameter 2 (7) 


Height of balloon in a zero pressure shape: 

Height = 0.748 • Diameter (8) 

Approximate average differential pressure due to buoyancy in a zero pressure design: 

\_ 

3 

(9) 


Af -0.834 
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Equations for a super-pressure pumpkin shape 


Diameter, top view, m: 


Diameter = 2.296 • 


0-75 M gas 

K n Pgas 


1 

J 
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Length of gore in a pumpkin balloon, m: 

L gore =1.312- Diameter (11) 

Surface area showing in a pumpkin shape: 

Aurf = 2.635 • Diameter 2 

or 2 (12a, 12b) 

A mrf =5.3 -Volume 3 


Height of balloon in a pumpkin shape: 


Height = 0.604 • Diameter 


Average differential pressure when in a super-pressure mode: 

i / n t 1 

_ _JE _ K z as ' 1 z as p 
Volume ^ 


(13) 


(14) 


If the unconstrained volume is greater than the maximum design volume Volume design , then the required average 
gas pressure is recalculated. If the balloon is a zero pressure design, then the gas will flow out the ducts with the 
elevated pressure according to well-established duct-flow equations. If the balloon is a pumpkin super-pressure 
design, the differential pressure can continue to rise until equilibrium is established or a vent is opened. 


P = P 

gas air 


+ A P 


(15) 


Illuminated Projected Area 

The illuminated projected area of a balloon varies with solar elevation angle ELV , and uses the top projected area 
A, op as the reference. Equation 16a was derived from information in Ref. 2. 

For a zero pressure shape: A proJecfed = A top * [0.9125 + 0.0875 cos(ti - 2 ‘ELV)] (16a) 

For a pumpkin shape: A projecled = A U)p * [0.8219 + 0.1781*cos(7t - 2 ELV)] (16b) 

For a spherical shape: A projected = A top ( 1 6c) 


B. Atmospheric Transmissivity 

The transmissivity of optical frequency irradiance thru the atmosphere follows a Beer’s Law format of 
exponential decay, the constants in the equation being derived experimentally for a clear atmosphere. A correction 
factor CF airmass is factored into the air mass ratio to account for fog and smoke, or for a different planet’s 
atmosphere. This is a modification to the original equation found in Ref.l. 

The Air mass ratio that the sun shines thru is calculated with the surface pressure P 0 : 1 


Air Mass = CF, 


f P N 

air 

\ P o J 


Vl229 + (614-sin (£LF)) 2 -614 -sin {ELV) 


(17) 


The transmissivity of a solar beam thru the atmosphere is then correlated with the air mass ratio: 


= 0.5.[ 


g-Q. 65- AtrMass ^- 0 . 95 -AirMas. 


„] 


(18) 


For the situation when the sun elevation angle is below the true horizon but above the apparent horizon, use this 
approximate formula to determine the air mass ratio (see Fig. 1): 


Air Mass 



1 + 


ELV 

Dip 


-70- 


ELV 

Dip 


( 19 ) 
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Figure 1. Dip of the horizon and air mass. 


I Sun - 


1358 


R 


AU 


C. Direct Solar Environment 

To generalize the solar environment for any planet, we must 
know the orbital eccentricity e and the mean orbital radius R AU . 
For small orbital eccentricities, one can approximate the true 
anomaly TA as being reasonably close to the value of the much 
easier to calculate mean anomaly MA: 

TA « MA = 2n — Day ^er (20) 

DaysPerYear 


The day number begins at perihelion (closest and hottest). 


Solar irradiance flux at the top of the atmosphere: 


’l + I. 


• cos(TA) 

2 

U-eJ 



Watts/m 2 


( 21 ) 


For Earth, R AU = 1, e = 0.016708, for Mars, R AU = 1.52368, e = 0.093400, and for Venus, R AU = 0.723, e = 0.0068 

At the balloon altitude Z, the direct solar irradiance is equal to the product of the top of the atmosphere value and the 
atmospheric transmittance: 

I SunZ = I Sun * T atm @2) 


D. Infrared environment 

Just as direct solar is attenuated for thickness of the atmosphere, the same is applied to ground and cloud IR that 
has to pass through a certain amount of atmosphere to reach the balloon. It is ASSUMED that the same coefficients 
in the equations from section C. apply to IR as well, even though there absorption lines for H 2 0 and C0 2 . 


Infrared diffuse flux at ground level with ground emissivity £ groU nd and ground temperature T grounc i (°K): 

I Rground ~ Aground * ® Aground Watts/m 


(23) 


As with visible light, the ground IR is attenuated with a similar formula to Eq. (18): 


W =1-716-0.5. 


- 0 . 65 -- 


- 0 . 95 - 


+ e 


(24) 


The ground-source diffuse IR flux at the balloon altitude Z: 

I RgroundZ ~ tflRground ^ atm I R Watts/m 2 (25) 


Deserts will have the greatest day/night temperature swings due to little water acting as heat capacity/storage. 
Somewhere on the order of 20°C swing is not unreasonable. Oceans should have the least day/night swing in 
surface temperature. In the software, BalloonAscent will use a user-defined maximum day and a minimum night 
temperature and sine wave between the two with a user-defined lag (in fractions of a day, 0.05 default). 

Here is a suggested list of ground emissivity : 

Desert; £ groU nd = 0.85 (data suggests great variation within a region, as low as 0.8) 

Average ground; £ groun d = 0.95 
SnOW, £ gr0 und ~ 0.98 

It is fortuitous that in the hottest places, the emissivity is low and in the coldest places the emissivity is high. 
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E. Albedo Environment 

There is a simple model relating surface albedo and 
albedo flux. The spacecraft EOS-TERRA has data on 
shortwave radiation flux average values and the 
calculated planetary albedo. The numbers are averaged 
monthly and so must include normal average cloud 
cover. The assumption here is that the albedo number is 
the total specular + diffuse solar reflection, but treated as 
diffuse. BalloonAscent uses the effect of solar elevation 
angle to achieve the instantaneous time varying albedo 
flux. One can see from Fig. 2 that the incoming beam of 
sunlight of flux strength Isun is diminished on its 
transmission to the ground thru the atmosphere. Since 
the solar angle is not directly overhead, the beam of light 
Figure 2. Albedo model at the top of the atmosphere [ s smeared out over a larger surface area, diminishing it 
further by the factor of sin(£Z,F). The measurements from spacecraft are from the top of the atmosphere, so this 
model under calculates for when the balloon is at low altitude. It is further assumed that atmospheric glow is part of 
the measured number. 

Albedo flux at the top of the atmosphere: 

Q albedo = Albedo ■ I Sun • sin (ELV) Watts/m 2 (26) 

F. Cloud Modifications 

The cloud fraction CF is defined as the fraction of the sky that is totally obscured by clouds. Note that clouds 
attenuate the direct solar when the balloon is under the cloud level. Planetary albedo is attenuated differently when 
below due to a one pass thru the clouds effect, and two passes through the clouds when the balloon is above the 
clouds. Earth IR is attenuated when above the cloud layer due to obscuration. When below the cloud level, there is 
no effect due to cloud albedo, but cloud albedo is counted only when above the cloud level. 

Solar beam radiation modified for clouds : 

If the altitude Z is below the clouds then direct solar flux q Sun = hunz * (1 - CF) watts/m 2 (27a, b) 

If the altitude Z is above the clouds then direct solar flux qsun = I Sun z 

Planetary diffuse IR radiation modified for clouds : 

If the altitude Z is below the clouds then qm P hnet = qiRgwmdz + CF • qmdoudz watts/m 2 (28a, b) 

If altitude Z above clouds then qiRpianei = qiRgmundz ‘ (1 - CF) + CF • qmdoudz 

Albedo modified for clouds : 

If the altitude Z is below the clouds then Albedo = Albedo Groun j • (1 - CF) watts/m 2 (29a, b) 

If altitude Z above clouds then Albedo = Albedo Groun( j • (1 - CF) 2 + Albedo C hud * CF 

G. Averaging the Film Optical Properties 

The optical surface properties are expressed as absorptivity a, transmissivity t, and emissivity £. The balloon 
surface is treated as an area-average between the properties of the balloon top (cap + shell) and the bottom (shell 
only), and then area-averaged again with the properties of the meridional load tapes. 

The view factor ViewFactor is the balloon surface area diffuse-radiant view factor of the planet. It is the ratio of 
balloon surface area that “sees” the planet surface divided by the total exposed balloon surface area. R p / a nei is the 
radius of the planet, meters, and is equal to 6371000 for earth. 

HalJCom^ = sin -1 (30) 

„ K Planet + Z _ 
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( 31 ) 


T/ . n , (l-cos (HalfCone le )) 

ViewFactor = - 

2 


H. Radiant Heat Loads on the Film 


The following equations will relate the radiant fluxes affecting the balloon skin, and the resulting absorbed 

radiant heat loads. As can be seen in Fig.3, the film is 
not only affected by the outer environment, but by the 
inner environment and the trapping of heat by multiple 
reflections. These multiple reflections raise the 
effective interior absorption. Recalling that 
reflectivity r = 1 - a - r, and the effective reflectivity 
becomes r effective = r + r I. 2 + r 3 + r 4 + r 5 +. . . 

From Kirchoff s law of radiation heat transfer, at 
any specific wavelength the absorptivity is equal to the 
emissivity while in thermal equilibrium. Even though 
this analysis is for transient conditions, we still follow 
this rule. This is important because absorptivity is 
generally given for solar short-wave frequencies, but 
emissivity is given for IR wavelengths. So if we are 
interested in the absorption of IR radiant energy, we 
use the given IR emissivity as our IR absorptivity ( a m 
= £). 

Absorbed direct sunlight heat: 

Qsun = a ■ ^projected ’ <lsu» ‘ [l + * ' 0 + r effec„j] Watts (32) 



Figure 3. Radiant flux sources affecting the balloon skin 


Absorbed albedo heat: 

Q Albedo = a • Aurf ' <1 Albedo ' ViewFactor • [l + r • (1 + r effectlve )\ Watts (33) 


Absorbed planetary IR heat: 


Qmptane, = <*,R ' Aurf ' VlRplane, ’ ViewFaCtOV • [l + T , R • (1 + 

Watts 

(34) 

Absorbed IR from the sky: 

Qmsky = aiR • Aurf • *1 1 Rsky * (l “ ViewFactor ) • [l + t ir ■ (1 + r effecli j] 

Watts 

(35) 

Absorbed IR self-glow from the interior: 

QlRfilm — O ' E ' (X J R * A sur j- • Tj}i m * (1 T y 'effective ) 

Watts 

(36) 

Emitted IR energy from both interior and exterior of the balloon skin: 

QiRout ® f • 2 • A sur f * T fllm 

Watts 

(37) 


I. Convective Heat Loads on the Film 

The convective heat transfer of the balloon involves transfers between the atmosphere and the exterior skin, and 
the lifting gas with the interior skin. The external convection is partitioned into free convection and forced 
convection. 

Air and Gas Properties 
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Dynamic viscosity of air 3 : 


Mai, = ■ 


1.458-10” 6 -T 


-6 rp \ .5 
air 


T n „ + 110.4 


Dynamic viscosity of helium gas 3 : 


=1.895-10 


-5 


( T \ 
gw 

v 273. 1 5 j 


0.647 


Conductivity of air 3 : 


k = 0.0241 


{ T \ 


V 273.15 j 


0.9 


Conductivity of helium gas 3 : 


kgas =0.144- 


T V ' 7 

1 gas 

^273.15 j 


Prandtl number for air, dimensionless 3 : 


iV, =0.804-3.25-10 -7\ 


Prandtl number for helium gas, dimensionless 3 : 


Pr gas =0.729- 1.6 -T gas 


N-s/m 2 (38) 


N-s/m 2 (39) 


Watts/m-°K (40) 


Watts/m-°K (41) 


(42) 

(43) 


External Free Convection 
Dimensionless Grashof number: 


Gr = 


Pair 'S' 


T -T 

1 film 1 ai , j 

T ■ u 2 

air r*air 


* Diameter 


Dimensionless Nusselt number for a sphere in free convection: 1 

Nu mr =2 + 0.45 -(Gr^-PrJ 25 


Free convection heat transfer coefficient for the external skin: 

Nu • k 
j-fr — 1 Y 1/1 atr 

free Diameter 


Watts/m 2 -°K 


(44) 


(45) 


(46) 


External Forced Convection 

When the balloon is ascending, the convection does not rely on differentially warm buoyant air for convection, 
but rather the vertical velocity forces the relative movement of air. First we calculate the Reynolds number, which is 
the dimensionless ratio of inertial forces to viscous forces; Re = Velocity ■ Diameter • p air / p air 

Forced convection heat transfer coefficient for the external skin: 3 

HC forceJ = ^ (2 + 0.41 • Re 0 55 ) Watts/m 2 -°K (47) 

Diameter ' 


Internal Free Convection 

Reference 3 has an equation for the heat transfer coefficient inside a sphere, which BalloonAscent uses unaltered, 
and is reproduced here: 
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internal ^ . 1 3 * k gag 


( _ 2 
Pgas 'g' 


T -T 

1 film 1 gas 


Pr. 


1 


gas 


rri 

* gas Pa 


Watts/m 2 -°K (48) 


Total Convective Heat Transfer 

For the external heat convection, the greater of the two values between free and forced convection is selected to 
equal HC exlernah The heat loads are thus: 


External convection heat load: 

QconvExt HC e xtemal ^ effective air 

Tfilm) 

Watts 

(49) 

Internal convection heat loss: 

Q Convlnt HC \ n t ema i ^effective film 


Watts 

(50) 


J. Lifting Gas Temperature Differential Equation 

The lifting gas rate of change of temperature is formulated on the adiabatic expansion response modified with 
internal convection interaction with the film. Whether from mass loss or purely a change of volume, or a 
combination of both, the gas will respond according to adiabatic expansion if the heat input Q Co nvim is zero. 


dT 


gas 


' Convlnt 


dt c„ ■ M, 


-+(r-i)- 


gas 


gas 


gas 


dp. 


gas 


dt 


°K/s 


(51) 


Expressed in terms of mass and volume change derivatives: 


dT 


gas _ 


' Convlnt 


dt M. 


+ (r-iK, 


gas 


dM 


gas 


1 


dt M 


gas 


dVolume 

~dt 


1 

Volume 


J 


°K/s 


(52) 


Where y= c/c v , and R gas / c„= y- 1 


K. Film Temperature Differential Equation 

The rate of change of the film temperature is derived from the simple transient energy-balance equation: 


dPfilm \Qsun " 1 " Q Albedo " 1 " GlRplanet " 1 " QlRsky QlRfilm QconvExt Qcomltu QlRoul ) 


dt 


C f ' ^ filmB 


°K/s (53) 


L. Valve and Duct Flow Differential Equations 

Flow thru a duct or valve depends on the differential pressure across the area interface, the cross sectional area, 
the density of the gas, and a discharge coefficient. The velocity thru such a duct or valve is: 

AP 

Velocity flow = C discharge - 2 ^ m/s (54) 

\ Pgas 


This leads to the differential equation: 

dM 


gas _ 


dt 


^ valve ^ discharge 


* ^d^valve Pgas 


kg/s 


(55) 


The mass flow can apply to either a duct flow situation as in a zero pressure balloon, or a valving situation as in a 
pressurized balloon. 
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M. 3-D Acceleration Differential Equations 

The total buoyant force produced by the balloon is the gross inflation, GI. With the average differential pressure and 
the gas temperature, the gas density is determined with the ideal gas law. The ambient air density is known, and 
with this information the specific buoyancy b is determined: 



^ S " ( Pair P gas ) 

N/m 3 

(56) 

The gross inflation: 





GI = b ■ Volume 

N 

(57) 


The aerodynamic forces are determined from the relative wind with respect to the balloon’s motion. If Ux, Uy, Uz 
are the balloon’s absolute velocity components, and Vwind x , Vwindy, Vwindz are the wind’s absolute velocity, then 
the following are the relative velocities: 

Vrelative x = Vwind x - U x 

Vrelative Y = Vwind Y - U Y (58a,b, c) 

Vrelative z = Vwind z - U z 

The magnitude of the relative wind: 

Vrelative = (Vrelative x 2 + Vrelative Y 2 + Vrelative z 2 y m/s (59) 

The magnitude of the drag force, where C d is the drag coefficient and the reference area is the top projected area: 

Drag = f • p mr • Vrelative 2 • C d ■ A, op N (60) 

The drag coefficient has in the literature for a sphere values ranging from 0. 1 to 0.47. However, a rising balloon is 
not a sphere, and more resembles an ice-cream cone which may have an equivalent drag coefficient based on the 
calculated spherical diameter of Q~ 0.8. The drag is now partitioned into its vector components: 

Drag x - Drag * Vrelative x /Vrelative 

Drag Y = Drag • Vrelative Y / Vrelative N (61a,b, c) 

Drag z = Drag • Vrelative z /Vrelative 


The differential equation governing the vertical motion of the balloon where M gross is the over-all mass of the 
balloon, less gas (suspended mass plus balloon material mass): 



GI + Drag z - M grms ■ g 

^ virtual 


m/s/s 


(62) 


The virtual mass is calculated with a virtual mass coefficient C virtuai . This is to take into account the mass of air that 
is necessarily dragged along with the balloon mass. Values assumed range from 0.25 < C virlua i < 0.5. 


M virtual = M gross + M g as + Cyirtua, * {Pair ’ Volume) kg (63) 


The horizontal accelerations: 



dt 


M virtual 



Drag? 

^ virtual 


(64a, b) 
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III. Software Interface 

The BalloonAscent interface is simple and intuitive. Figure 4 shows the main form in the background with a top 
row menu bar, which includes some handy balloon flight calculators. The sandwiched form is the data page, and it 
is subdivided into 13 categories, each accessible by clicking on the appropriately labeled command button. The data 
fields are made visible, sometimes with a new pop-up form, as is shown in the ambient environment data card. On 
the top row menu bar there are location templates to load coordinates and dates from the typical NSBF launch sites, 
as well as templates for typical NSBF zero pressure balloons. The ambient environment is defined such that it can 
vary four different times during the simulation. The maximum daytime and minimum nighttime temperature for a 
given area is defined, with a sine-wave variation between the two implemented automatically. The atmosphere can 
also be modified with a diurnal delta temperature variation and lag time, which will effectively compress the 
atmosphere slightly at night. The ground emissivity and albedo is defined, as well as the cloud albedo and cloud 
fraction, which just about sums up the typical float environment. As a bonus, which also helps to model off-world 
situations, the radiant flux environment (solar, albedo, IR up-welling, 1R down-welling) can be directly defined 
using Beer’s law exponential decay coefficients. The software has a calculator on the ambient environment data 
page to calculate the coefficients from planetary data that you might happen to have. 


JOJJL] 


f Rodgers Balloon Ascent and Thermal Simulator 


dick for Data Page Print TItIs Page HELP CALCULATORS EXIT PROGRAM 

Press button to graph variable 
_J Bulk Film Temperature 

I Bulk Gas Temperature 

_J Air Temperature 
I Differential Pressure. Apex 


Data file = C:\DATA FILESW/isualBasic 5 Files BalloonAscent l\Ascent Data Files \6mcf\6 mcf Ft Sumner AugSOBDA 


CONDITIONS AT THE SIMULATION END 
Total Ballast Dropped, kg 657 

Gas Mass, kg 377.1 
Average Gas T emperature, K. 218 

Average Film Temperature, K 218.3 

Air Temperature, K 233.9 

'/ertical Speed, i 
Altitude, m 

Latitude. Longitu Read from Fie Save to Fie Load Balloon Templates Load Locations 

Distance travele 
Sun Elevation, d 
Apex Differential 
Launch Stress In 
Launch Free Lift 
Required Cop Li 
Accumulated Te 


P * 


Balloon 

Characteristics | 

Launch Location 
and Time 

Ballast Control 

Optical Properties 

Winds Aloft 
f(altitude) 

Lifting Gas 

Planetary Data j 

Material Properties 

Ambient 

Environment 

Valve Control 

Vent Ducts 

Atmosphere 

Miscellaneous 

Settings 


dai i rtnu ruADArrcDicnrc 





Heat Loads into 
Qsun 0 
Qolbedo 0 
QIR planet 2f 
QIRsky 527 
QIR 'film glow 1 1 
QconvAir 23 
QconvHelium 

Radiant Source ' 
Average qsun fli 
Average qalbed 
Average qlFtplai 
Average qIRsky 
Average Albedc 
Source qlRgrour 
Source qIRclouc 
Ground IRAtteni 

Average skin ab 
Average skin IR 
Average skin so 

Convective Hea 
HC internal, watt: 
HC external, watt 


[445 

jl363 


Mass of Gas. kg 


* AMBIENT ENVIRONMENT 


|1143 

f?77448 


1^2 

Max Day Temperature of 

Local Geography, K [318 [285 


[22 


[300 


Min Night Temperature of f^gg f^jjg |298 
Local Geography. K 


|6mcf pumpkii 

BALLO- 


Emissivity of Local 
Geography 


Ground-Only Albedo 


[~ Press to 


Write Output to T ext File 


Click to Start Simulation 


Cloud-Only Albedo 


Cloud Fraction, 
obscured sky 


|0 

[oT 


| ! 

|o 05 


|12 

|0.95 


|12 

fd.23 


J12 

ET 


J 12 

loir 


Right Time, hrs 
Day Temp, K 
Night Temp.K 

Right Time, hrs 
Temperature. K 

Flight Time, hrs 


|22 

foT" 
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Figure 4. Program 1 st form, main data page, and ambient environment data card 
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IV, Example Results 


A. Adiabatic Test Runs 

To test whether or not the simulation was correctly determining the gas temperature during ascent or during a 
period of venting, the internal convective heat input to the gas Qconvim was set to zero. One would then expect the 
gas to behave in an adiabatic fashion according to the following equation: 

/ \r- 1 


t 2 =t, 


Pi 


(56) 


\Px) 


Several runs were set up, the first an ascent with no venting, the second a pressurized pumpkin venting off 
pressure at a constant altitude. In both cases BalloonAscent followed the adiabatic expansion Eq. (56) exactly. 

Considering that helium is a mono-atomic gas (no rotational or vibratory degrees of freedom) with weak Van de 
Walls forces when warm, it should follow the ideal gas laws assumed here quite well. Expansion thru small nozzles, 
where Thompson-Joule effects will deviate temperatures greatly from adiabatic assumptions, is not credibly 
applicable to ballooning except in situations during filling from helium bottles. This software does not address that 
portion of balloon operations, nor is it necessary for determining ascent or float behavior. 

B. Zero Pressure Balloon Simulation 

The 40 “heavy” zero-pressure balloon example was simulated here for a flight from McMurdo station Antarctica. 
Such balloons have seen diurnal variations in altitude on the order of 4500 meters, where as the simulation shows a 
larger variation on the order of 8000 meters (see Fig. 5). The change in altitude is caused by the variance in albedo 
flux, as can be seen in Fig. 6. The most logical explanation is that snow and ice have specular as well as diffuse 
albedo components, which will add complications to the simple diffuse albedo model of Eq. (26). Future versions 
will therefore include a specular albedo flux as a function of solar elevation angle and zenith surface reflectivity. 

Figure 7 shows the gas, film, and air temperatures, and Fig. 8 is a ground track on an orthographic projection 
centered over the average location. 


jJOlxj 
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pOmd heavy Polw Right 

Figure 5. Zero pressure balloon altitude simulation 
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Figure 6. Zero pressure balloon environmental flux 
levels 
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* Graph Form 


i jnu 
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Figure 7. Zero pressure balloon temperatures 



Figure 8. Zero pressure balloon ground track 


C. Super-Pressure Balloon Simulations 

The super-pressure example simulation is a 13 million cubic foot pumpkin with 245 gores, launched from Alice 
Springs Australia. It simulates the hot day case of traveling over the desert, and a cold case that night over a cold 
storm in the Pacific Ocean followed the next day by a typical mid-latitude warm ocean day. Figure 9 shows the 
variation in altitude; indeed the balloon looses all super pressure during the cold storm night and levels out at a 
lower altitude. Ballast could have been dropped, but to achieve 100 days of flight, an occasional altitude drop if not 
too severe is considered acceptable. 
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Figure 9. Super pressure balloon altitude variation 


-!□! > 



Elapsed Time, hr 


48 


Figure 10. Super pressure balloon temperature 
variation 
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Figure 11. Super pressure balloon ascent velocity Figure 12. Super pressure balloon differential pressure 


D. Mars Balloon Simulation 

The two things to understand about Mars balloons are that scale-wise they fly at the bottom of the atmosphere. 
The atmosphere is thin, so the bottom offers the most buoyancy. This also forces the design payload to be small 
compared to the mass of the balloon material. The other consideration is the large diurnal temperature variation of 
the surface and atmosphere, reportedly 30 °K. This has the effect of compressing down the atmosphere at night, 
which leads to a very curious phenomenon regarding balloon flight on Mars if using a super-pressure design. The 
idea of the pressurized balloon is that so long as some internal pressure is maintained, the displaced volume is a 
constant. This means that the system density (total mass including gas + total displaced volume) is constant, and as 
all balloons, will match the atmospheric density when in float equilibrium. A Mars super-pressure balloon will find 
itself at night lighter than the surrounding air that has compressed down, and as a result will float upwards at night, 
and sink back down during the day, as shown in Fig. 13. Figure 14 shows the resulting radiant flux environment at 
the balloon’s altitude. 
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Figure 13. 575 cubic meter Mars balloon altitude 
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Figure 14. 575 cubic meter Mars balloon environment 
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V. Conclusion 

BalloonAscent is on its 32 nd version and still evolving. It is extremely user-friendly with intuitive data entry 
interfaces. The environment is defined with measurable parameters, and transmission models are used to propagate 
the environment to the balloon’s altitude. Envisioned improvements in the future are to have built-in sigma shapes 
such that stub-ducts can be correctly modeled. A stub-duct is a short duct that is opened on command with a squib 
in order to lose a specific amount of gas such that a specific amount of altitude is lost. The other improvement has 
been mentioned already, to improve the albedo model for polar flights by including a specular component to the 
model (currently only diffuse is assumed). Another improvement is to consider a simple down-welling sky IR 
model as a function of altitude, which could add as much as 160 watts / m 2 at the launch altitude. This effect is 
admittedly very small, and at float the effect is nil. 
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